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ABSTRACT 

We present the results of the Velocity Coordinate Spectrum (VCS) technique to calculate the velocity 
power spectrum of turbulence in the Small Magellanic Cloud (SMC) in 21cm emission. We have 
obtained a velocity spectral index of -3.85 and an injection scale of 2.3 kpc. The spectral index is 
steeper than the Kolmogorov index which is expected for shock-dominated turbulence which is in 
agreement with past works on the SMC gas dynamics. The injection scale of 2.3 kpc suggests that 
tidal interactions with the Large Magellanic Cloud are the dominate driver of turbulence in this dwarf 
galaxy. This implies turbulence maybe driven by multiple mechanisms in galaxies in addition to 
supernova injection and that galaxy-galaxy interactions may play an important role. 

Subject headings: ISM: structure — MHD — turbulence, SMC 


1. INTRODUCTION 

Magnetohydrodynamic (MHD) turbulence is ubiqui¬ 
tous in the interstellar medium (ISM) of galaxies. It is an 
essential element in many physical processes, such as the 
formation of stars (Federrath & Klessen 2012; Burkhart, 
Collins & Lazarian 2015) and transport of cosmic rays 
(Lazarian & Yan 2014), and is a large fraction of the en¬ 
ergy budget of the Milky Way (see Elmegreen & Scalo 
2004). The evidence for the existence of turbulence in 
the ISM is overwhelming and includes the “Big Power 
Law” and “Extended Power Law” of the electron den¬ 
sity fluctuations in Armstrong et al. (1994) and Chep¬ 
urnov & Lazarian (2010), respectively, fractal structure 
of molecular clouds (Stutzki et al. 1998), and intensity 
fluctuations and broadened linewidths in channel maps 
(see Crovisier & Dickey 1983, Green 1993, Stanimirovic 
et al. 1999, Deshpande, Dwarakanath & Goss 2000). De¬ 
spite the importance of ISM turbulence for a wide range 
of physical processes in galaxies, little is known about 
it due to the difficulty of measuring the velocity power 
spectrum of turbulence which provides information on 
the injection and dissipation scales, energy amplitude, 
and spectral slope. 

For studies of turbulence, velocity fluctuations re¬ 
flected in position-postion-velocity (PPV) radio cubes 
are critical as they trace the underlying fluid behavior. 
Along with density fluctuations they contribute to the 
intensity fluctuations in PPV space. As such, many tech¬ 
niques for analyzing PPV fluctuations in the context of 
turbulence have been proposed including velocity cen¬ 
troids (Esquivel & Lazarian 2005a,b; Ossenkopf et al. 
2006; Esquivel et al. 2007), hierarchical tree methods 
(Goodman et al. 2009; Burkhart et al. 2013), principal 
component analysis (Heyer & Schloerb 1997) and the 
A variance technique (Ossenkopf et al. 2008a,b) How¬ 
ever, these techniques face difficulty to reliably separate 
out the underlying velocity and density contributions in 
PPV space in the supersonic limit. This is essential for 
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obtaining the velocity power spectrum of turbulence in 
astrophysical settings. 

For the purpose of obtaining the velocity power spec¬ 
trum the Velocity Channel Analysis (VCA, Lazarian & 
Pogosyan 2000, 2004, henceforth LP00, LP04) and Veloc¬ 
ity Coordinate Spectrum (VCS, see Lazarian & Pogosyan 
2000, 2006, 2008, henceforth, LP00, LP06, LP08) have 
been developed based on analytic expression for the be¬ 
havior of the correlation function (which is isomorphic to 
the power spectrum) in PPV space. A concise review of 
these extensive works is given in Lazarian (2009). 

The VCS and VCA techniques are entirely data driven 
(i.e. comparisons or scaling to numerical simulations are 
not required) and are complimentary to each other as the 
success of the VCA application depends primarily on the 
spatial resolution while the VCS depends primarily on 
the spectral resolution. This is because the VCA tech¬ 
nique is concerned with taking the spatial power spec¬ 
trum of the two-dimensional slices of PPV while the VCS 
studies the fluctuations of PPV intensities along the ve¬ 
locity coordinate of the PPV cube. Both techniques have 
been applied to multi-wavelength observations as well 
as undergone numerical tests (Stanimirovic & Lazarian 
2001; Chepurnov et al. 2010; Chepurnov et al. 2009) 

In this paper we apply the VCS to the Small Magellanic 
Cloud (SMC) in 21cm emission. The SMC is a dwarf ir¬ 
regular galaxy in the Local Group which has a gas rich 
ISM environment and is metal poor (see, Stanimirovic 
et al. 1999, henceforth known as SX99). The SMC is 
nearby (60 kpc; Westerlund 1991) and is close enough 
to have resolved scales HI from ss 4.5kpc down to 30pc. 
This is in contrast to clouds in the Milky Way where dis¬ 
tance determination is relatively uncertain and special 
geometrical considerations must be taken into account 
for application of the VCS as was done in the applica¬ 
tion of the VCS in high latitude HI data (Chepurnov et 
al. 2010). Another advantage of investigating the SMC 
is that it has been used as a test bed for numerous turbu¬ 
lence studies. HI observations of the SMC, obtained us¬ 
ing the Australia Telescope Compact Array (ATCA) and 
the Parkes telescope (SX99), have already been studied 
with the VCA technique (Stanimirovic & Lazarian 2001) 
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as well as the traditional power spectrum. Other statis¬ 
tics were applied to the SMC by Burkhart et al. (2010) 
which provide information on compressibility, (i.e. sonic 
Mach number) such as the PDFs and spatial power spec¬ 
trum. In addition, the topology of the SMC and its rela¬ 
tion to the turbulence environment was quantified using 
the genus statistic in Chepurnov et al. (2008). Because 
the SX99 SMC data set is so well studied, it is a perfect 
candidate to test and compare the newer VCS technique. 
Thus by studying this data set we can acquire new infor¬ 
mation, test and confirm past results, and validate the 
promise of these tools for further use in other observa¬ 
tional studies. 

This paper is organized as follows: in Section [2] we de¬ 
scribe the SMC HI data set used, in Section [3] we review 
the VCS technique, in Section 2] we detail the results of 
our analysis and finally in Sections and [5] and 2] we discuss 
the implications of our findings and draw our conclusions. 

2. THE SMC IN HI: TEST GALAXY FOR TURBULENCE 
STUDIES 

The SMC HI data set used in our analysis combines 
both single dish (Parkes, Stanimirovic et al. 1999) 
and interferometer (ATCA, Staveley-Smith et al. 1997) 
Both sets of observations were then combined (see Sta¬ 
nimirovic et al. 1999), resulting in the final HI data 
cube with angular resolution of 98", velocity resolution 
of 1.65 km s -1 , and 1-a brightness temperature sensi¬ 
tivity of 1.3 K. The data has a continuous range of spa¬ 
tial scales between 30 pc and 4 kpc. The velocity range 
of the observations is 90-215 km s -1 . For more details 
about the ATCA and Parkes observations, data process¬ 
ing, and data combination (short spacings correction) 
see Staveley-Smith et al. (1997) and Stanimirovic et al. 
1999. HI self-absorption we corrected for following the 
scheme of Stanimirovic et al. 1999. The correspondent 
HI intensity image is shown in Figure [I] 

3. OVERVIEW OF THE VCS TECHNIQUE 

The suggestion to analyze the data calculating the 
spectrum along the velocity axis of PPV cubes was first 
made in LP00 and elaborated in LV06. It was only later 
that the technique was applied to actual observational 
data (see Padoan et al. 2009, Chepurnov et a. 2010). 

LP06 mostly dealt with the asymptotic scaling of tur¬ 
bulence and the observational study by Padoan et al. 
2009 is an illustration of the application of the VCS to 
molecular clouds. The VCS was applied to galactic HI 
data in Chepurnov et al. (2010) not only to obtain the 
spectral slope, but also the injection scale of the turbu¬ 
lence as well as other pieces of useful information. Here 
we will review the main points of the VCS technique fol¬ 
lowing the way as it is presented in Chepurnov & Lazar- 
ian (2009). 

Below we briefly present the idea of the VCS technique. 
A more detailed discussion, including the analysis of the 
effects of absorption is given in Lazarian & Pogosyan 
(2006). The contributions to the PPV power spectrum 
is are coming the power spectrum of velocity, P v oc k ~ av , 
and the power spectrum of emissivity, P e oc k ~ ae , where 
k is the wavevector and k oc l/l where l is the spatial 
scale Power spectra P v and P e , respectively, determine 
the turbulent cascade of energy and the distribution of 
density in the interstellar medium. The VCS technique 


obtains the underlying velocity power spectrum by calcu¬ 
lating the ID Fourier power spectrum of intensity fluctu¬ 
ations (i.e. Pi ( k v )) along the velocity axis of a PPV cube. 
The VCS is based on the theory that relates the fluctu¬ 
ations of intensity in the PPV scale to the underlying 
fluctuations of velocity and density. The procedure that 
we employ is based on the fact that the Pi(k v ) depends 
on the angular resolution of the telescope beam (see Fig¬ 
ure 1). By investigating how Pi changes when changing 
from high to low angular resolution, and by utilizing the 
information from intensity map (or column density map) 
we can evaluate the velocity power spectrum. 

The spectral line signal that we deal with is the HI 
intensity measured at velocity i>o and at a given beam 
position e: 

S(e,v 0 ) = J w(e, r) dr e(r)f(v r (r) + v r r eg (r) — u 0 ), (1) 

where e is normalized emissivity,/ denotes the convolu¬ 
tion between the spectrometer channel sensitivity func¬ 
tion and the temperature dependent Maxwellian distri¬ 
bution of velocities of emitting particles. In addition, 
v^ eg is a line-of-sight component of the regular galactic 
velocity shift, v r is the line-of-sight component of the 
random turbulent velocity. The window function w is 
defined as: 

w(e, r) = -^(e, r)w e (r), (2) 

where Wb is an instrument beam, pointed at the direction 
e, which depends on angular coordinates r, while w E is 
a window function defining the extent of the observed 
object. 

The Fourier transform of a spectral line is 

S(e, k v ) = ± S(v 0 )e~ ik ^° dv 0 

= f{k v ) J w(e,r)dr • ( 3 ) 

e(r) exp(—?'fc„(ty(r) + v r r eg (r))). 

If we correlate S in two directions, pointed by ei and 
e 2 , we get the following measure 

A(ei, e 2 , k v ) = (s(ei,k v )S*(e 2 ,k v )^ 

= f 2 (k v ) J w(e i,r) dr f w(e 2 , r') dr' ■ 
(e(r)e(r')) (exp(-ik v (v r (r) - iv(r')))) • 
exp(—ik v (v^ eg (r ) - V s (r'))), 

( 4 ) 

where (...) denotes averaging and where we assumed that 
gas velocity and emissivity are uncorrelated. 

Equation 2] gives us the emissivity correlation function: 

Ce(r-r') = (e(r)£(r')). (5) 

We provide averaging assuming that the velocity statis¬ 
tics are Gaussian^ 

(exp(-ik v (v r (r) - v r ,(r')))) 

= ex P (“%■ (K(r) - iv(r')) 2 )) • 

We assume that the beam separation and the beam width 
are both small enough that we can neglect the difference 

3 We assume that the velocity field is Gaussian. The latter is 
true to high accuracy (see LP00). 


Velocity Coordinate Spectrum of the SMC 


3 



Fig. 1.— SMC integral intensity map, with the contour surrounding region with the signal level being more than 20% from the maximal 
intensity, selected for the VCS computation. The spatial scaling is in parsecs. 


between v r and v z (we consider z-axis to be a bisector of 
the angle between beams). We also assume that v r z eg (r) 
depends only on z : 

v r / g (z) = b(z-zo) + v™ g , (7) 

where b characterizes regular velocity shear. The case of 
velocity shear b arising from Galaxy rotation is consid¬ 
ered in detail in LP00. 

Let us introduce the velocity structure function: 

D vz { r - r') = (K(r) - u z (r')) 2 ) , (8) 

then Eq. 0 can rewritten as: 

K(ei,e 2 ,k v ) =f 2 (k v )- 

f w(e i, r) dr J w(e 2 , r') dr' C £ (r — r')- 

exp (-^D vz {r - r') - ik v b(z - z')j . 

( 9 ) 

Further transformations result in 
K(e i,e 2 ,k v ) = f 2 (k v ) J g(ei, e 2 , r) dr ■ 

C e (r) exp (~^fD vz { r) - ik v bz^j , 

where geometric factor g is defined 

g(ei,e 2 ,r)= w(ei,r')w(e 2 ,r' + r) dr', (11) 


It is easy to see that the measure K given by Eq. m 
depends both on the angular distance between the vec¬ 
tors ei and e 2 and the velocity wavevector k v . If we 
integrate over k v we are left with the 2D power spec¬ 
trum that is the basis of the VCA technique. On the 
other hand, if we take coincident vectors ei and e 2 to 
obtain a one-dimensional spectrum P 1; we arrive at the 
VCS technique that we employ in this paper. 


Pi{k v ) = K(e u e 2 ,k v )\ &1=&2 = f 2 (k v ) f g(r) dr ■ 

C e ( r) exp (y-^fD vz { r) - ik v bz^j . 

where g(r) = g(e, e, r) and e is a beam direction. 

Eq. (fl2l) has two asymptotic spectral regimes that de¬ 
pend on beam size (see LP00 and CL06) and on whether 
the density spectral slope is shallower than -3 or steeper 
than -3. The gist of our approach is that we do not use 
the asymptotic regimes but rather calculate PI directly. 

First of all, we derive the window function. SMC is a 
remote emitting structure and therefore the correspond¬ 
ing window function w e is non-zero only in some vicin¬ 
ity of a distant enough point zq- We also assume, that 
picture-plane extent of the object is larger than the beam 
size, so the R-dependence of w e can be ignored. Then 
we obtain 


w;(r) 


~ 2 W b (9)w e (z ) 


and / is a Fourier transform of the channel sensitivity 
function /. 


(13) 
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Eddie modes: 

1 - low resolution 

2 - high resolution 

3 - intermediate 


Instrument beam 



Emitting structure 



Fig. 2.— The relation between the VCS technique and the effective beam resolution. Eddies that lie resolved within the telescope beam 
are labeled 1 and are considered part of the “low-resolution” mode, eddies with size larger than the telescope beam are labeled 2 and are 
considered part of the “high resolution mode”. 


If we introduce 


presented as a sum of variable and constant parts: 


W b (R) = \w b (-), f W b (R ) dR = 1 (14) 

z o \z 0 J 

window function becomes 


w (r) = W b (R)w e (z) (15) 

and the geometric factor g( r) is the auto-correlation of 
it. In our calculations we use Gaussian approximation 
for both W b (R) and w e (z). 

Our goal is to express D vz through the velocity spec¬ 
trum. If we assume that v(r) is primarily solenoidal with 
power-law power spectrum having cutoff at large scales, 
the velocity powe r spectr u m ca n be written in as follows 
(see, for example. iLesieuHflQQll l: 


Ei 


= it*"' 


hkj \ 
k 2 ) ’ 


(16) 


where E 0 2 is the velocity power spectrum amplitude, a v 
is the velocity spectral index, k 0 = 2n/L v is the cutoff 
wavevector bound with the injection scale L v , and i and j 
are the component indexes. Then D vz can be represented 
as 

D vz (r) = 2 J dk (1 - e ikr )ziZjFij(k), (17) 


C e ( r) = J A\Fl{ k)e ikr dk + el (18) 

where A 2 .F 2 (k) is the emissivity power spectrum and £o 
is the mean emissivity. Here A 2 is the spectrum ampli¬ 
tude. 

For our modeling we need the value of 6/ = A 2 /e§. 
After some algebra we can derive for it the following ex¬ 
pression: 

8 *= (^?) |W ? 6,o(K)| 2 |ui £ ,ofe)| 2 F 2 (k)dk) 

where AS/So is the relative deviation plane of sky in¬ 
tensity, that is equal to 0.402 for the SMC intensity map 

within the selected region. W b ,o and w £ , o are the window 
function Fourier transforms normalized to their values at 
zero wavevector. Then we can evaluate the emissivity 
correlation function as follows: 

C e (r) ~ J S/F/{k)e lkr dk +1 (20) 


where 


^ e 2 (k) = 


_ 0 


( 21 ) 


where the summation over repeating indexes is assumed. 
We point out that Eq. [16] contains the injection scale 
L v = 2ir/ko and the velocity power spectrum amplitude 
Eq 2 . Thus the injection energy and injection scale can also 
be found from the VCS in addition the power spectral 
index. 

In addition, we need to estimate the emissivity cor¬ 
relation function for our modeling. This part can be 


For our calculations we assume that ko is the same as 
in Eq. [16] In what follows, we use the mathematical 
formalism above to analyze the observational data. 

4. RESULTS 

The two-dimensional spatial power spectrum of the HI 
column density image was performed by Stanimirovic & 
Lazarian (2001). They found that a e oc 3.3. Using a 
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distance to the SMC of 60kpc and thickness (size) of 
the SMC 4kpc (SX99) we apply the VCS procedure as 
outlined in the previous section. We first calculate the 
power spectrum over the velocity coordinate P\ for three 
resolutions: 0°.025, 0°.05 and 0°.l. The resulting spectra 
are shown in Figure [3] We note that our analysis is 
mostly biased towards the cold gas in the SMC and that 
the spectrum measured here is generally of the cold gas, 
which tends to have higher sonic Mach numbers than 
the warm HI in the SMC (Burkhart et al. 2010). The 
fitting of the spectrum, which is overplotted in Figure 
H is performed with Mathematica FindMinimum, which 
implements the algorithm of steepest descent. 

Applying the VCS to the SMC data we report the fol¬ 
lowing parameters: 

• a v = -3.85 ±0.02 

• a e = —3.3 (Stanimirovic & Lazarian 2001) 

• L v = 2.3lg;4 kpc 

• Vt.urb = 6.7_o.4 kill/S 

• M cold = 

The mach number for the cold phase M co id is calcu¬ 
lated assuming the cold phase temperature T co i d = 173K. 
We note that the P\ fitting did not include the first 3 
points, as they are affected by the regular motions (see 
below). 


4.1. SMC Regular Motions and the VCS 

An important consideration for the analysis of veloc¬ 
ity statistics in a galaxy such as the SMC is the con¬ 
tribution to the turbulence power spectrum from regu¬ 
lar motions. Regular motions in the SMC were stud¬ 
ied in Stanimirovic, Staveley-Smith & Jones (2004) who 
found a complicated velocity field resulting from several 
expanding HI shells in addition to the galactic rotational 
velocity. 

Let us estimate the possible impact of these motions. 
The velocity term of the VCS for arbitrary regular ve¬ 
locity is: 

exp(-|U(r - r') - ik v (v^( r) - < es (r'))) (22) 


Comparing the arguments for the turbulent and regular 
velocities and accounting for the fact that the spatial 
dependency of D z is three-dimensional, while that of v r z e9 
is one-dimensional (due to the relatively narrow beam) 
we can estimate the boundary k bound which separates 
regular motion domination region from the region where 
turbulent motions dominate as follows: 


k 


bound 

v 


/Av r z eg 

V «t„ r 6 


(23) 


We estimate An” 6 ® as the separation between large scale 
regular components of the average line-of-sight velocity. 
According to our calculations, Ad 69 = 48000 m/s and 
Vturb = 6700 m/s. Substituting this to the previous ex¬ 
pression we finally have: 


\ogk b v ound ss -3.6 (24) 


We note that this estimation fits well to the approximate 
scale at which our VCS fitting deviates from the data. 

5. DISCUSSION AND CONCLUSIONS 

The velocity power spectrum is a critical parameter of 
turbulence to obtain from ISM observations as it gives 
true dynamical insight into the process of turbulence. 
This is not the first time the velocity power spectrum 
has been estimated for the SMC galaxy as Stanimirovic 
& Lazarian (2001) used VCA to recover the velocity 
slope. While our estimate for the density spectrum in 
this paper agrees well with the value obtained using the 
VCA technique in Stanimirovic & Lazarian (2001) our 
estimate of the velocity spectrum is substantially differ¬ 
ent. Indeed, Stanimirovic & Lazarian (2001) found the 
three-dimensional velocity power spectral slope to be - 
3.4, which, in fact, is an unusual value for MHD turbu¬ 
lence (see Brandenburg & Lazarian 2013 for a review). 
In contrast, we find a much steeper slope for the velocity 
power spectrum of -3.85. This slope is predicted from 
numerical simulations of high Mach number turbulence 
(e.g. Kritsuk 2007; Kowal & Lazarian 2010; Collins et al. 
2012) as well as the sonic Mach number analysis of the 
SMC done in Burkhart et al. (2010). The steep veloc¬ 
ity slope in a compressible medium also harkens to the 
predictions of a steep velocity slope given in the Burgers 
shock dominated turbulence model. 

We believe that fitting of integral expressions that 
we employed here provides a more accurate approach 
compared to the asymptotic approach applied in Sta¬ 
nimirovic & Lazarian (2001). The observed spectrum 
of fluctuations is limited by the resolution of the data 
and, given the limits of sensitivity, the uncertainty of the 
spectral index evaluation may be substantial. For in¬ 
stance, if we consider measuring the slope of the spectra 
of 2D channel maps within VCA, we should be aware 
for a given channel thickness the slices may show “thick” 
slice asymptotics for smallest eddies, while for large ed¬ 
dies the slices will exhibit “thin” asymptotics. Thus the 
spectrum at small wavenumber may have a spectral slope 
which differs from that of large wavenumbers and fitting 
a single power law may be misleading. A possible ap¬ 
proach to dealing with this issue is choosing the channel 
thickness of the order of the thermal line width. This 
however may require much higher signal-to-noise data. 
As a result, a similar procedure of fitting the spectra by 
using integral expressions and varying the thickness of 
the channel maps may be useful for the VCA practical 
analysis of the data. We believe that this should reconcile 
the data of the VCA with our results. Nevertheless, we 
want to point out that the compressible and incompress¬ 
ible components of velocity field contribute differently 
into the measures available to the VCA and VCS. This 
may be another reason for the disagreement that we face. 
The issue requires further studies. 

What drives the turbulence in the SMC? The fitted 
value of the 2.3 kpc injection scale suggests two domi¬ 
nate possibilities. One possibility is that turbulence is 
stirred by a large number of expanding HI shells found 
in the SMC. The shell sizes range from «30 pc to «2 
kpc. Hence, one scenario could be that the largest shells 
drive the turbulent cascade down to the smallest ob¬ 
served scales or that the shells fragment and drive in¬ 
stabilities. Turbulence may also be driven in the SMC 
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Fig. 3.— Fitting of the model VCS spectra to P\ of the SMC for different effective resolutions. We do not include in our fitting the first 
three points as they are affected by the regular motions. 


by the last SMC-LMC encounter. This would agree with 
the predictions made by simulations presented in Besla et 
al. (2012) and also with the analysis of Goldman (2000). 
We did not find evidence for a dissipation scale present 
in the HI data however the scales resolved in our data 
only proceed down to 30pc and thus the dissipation scale 
in the SMC may be smaller than this scale. This agrees 
with the line width size relation of the SMC, which as 
been observed down to scales smaller than 30pc (Rubio, 
Lequeux & Boulanger 1993). 

In light of the VCS analysis we can estimate the en¬ 
ergy density of turbulence in the SMC. Given a turbulent 
velocity at the injection scale of v tur b=6700 m/s and an 
average ambient HI density of 0.3-1.9 atoms/cm 3 (SX99), 
then the energy density of turbulence in the SMC would 
range between w tU rn = 0.07 — 0.44 eV/cm 3 , as com¬ 
pared to « lev/cm 3 for the Milky Way. As a compar¬ 
ison, the energy density in cosmic ray electrons in the 
SMC has been estimated by the Fermi telescope to be 
=0-0036 eV/cm 3 , which is about 60% of the local 
cosmic ray energy density in the Milky Way (Abdo et al. 
2010), and the energy density in the cosmic ray protons is 
estimated to be wcr p =0-15 eV/cm 3 (Abdo et al. 2010). 
Mao et al. 2008 estimated a total 3D field strength in 
the SMC (assuming equipartition with cosmic rays) of 
1.7±0.4/rG which suggests a magnetic energy density of 
«0.12±0.06 ev/cm 3 . The Milky Way appears to have 
an equipartition of energy between turbulence, magnetic 
fields, and cosmic rays with each component carrying 
about 1 ev/cm 3 . The SMC is a dwarf galaxy (with dy¬ 
namical mass of 2.4xlO 9 M@) and thus it is unsurprising 
that the overall energy density in turbulence, magnetic 
fields, and cosmic rays, is significantly less. 


The VCS technique, which as been numerically tested 
in Chepurnov et al. 2009, has a bright future ahead 
in terms of measuring the properties of turbulence in a 
number of interesting environments. The technique thus 
far has been applied primarily in the context of 21 cm 
data but there is no reason not to apply it to the wealth of 
molecular line tracers such as CO as was done in Padoan 
et al. (2009). Accurate measurements of the turbulence 
energy budget of molecular clouds and the sonic Mach 
number will have important consequences for theories 
of star formation (i.e. see Krumholz & McKee 2005; 
Padoan & Nordlund 2002; Hennebelle & Chabrier 2008; 
Federrath & Klessen 2012, and references therein). The 
VCS does not need to be limited to ISM or galaxy studies 
but also can be applied to spectral lines as measured 
in galactic halo environments or galaxy clusters using 
absorption lines or x-ray emission lines with high enough 
spectral resolution. 

Our results in this paper are in agreement with 
the study of Mach numbers of turbulence in SMC in 
Burkhart et al. (2010). This illustrates well the syn¬ 
ergy of using different techniques to investigate the pa¬ 
rameters of turbulence. A new technique of studying 
turbulence spectra of magnetic fluctuations using syn¬ 
chrotron polarization fluctuations was recently developed 
in Lazarian & Pogosyan (2015). This technique makes 
use of a statistical description of synchrotron statistics 
and applicable to describing fluctuations for arbitrary 
spectral index of cosmic rays. It will be advantageous 
to study and compare the statistics of velocity and mag¬ 
netic field in order to obtain a comprehensive picture of 
turbulence in SMC and other galaxies. 


6. CONCLUSIONS 







Velocity Coordinate Spectrum of the SMC 


7 


We apply the Velocity Coordinate Spectrum technique 
to obtain the velocity power spectrum in the SMC galax¬ 
ies in 21cm emission. We find a steep power spectral 
slope for velocity which is in agreement with the expec¬ 
tations for supersonic turbulence. In particular, we ob¬ 
tain the following properties of the turbulence velocity 
spectrum in the diffuse gas of the SMC: 


• The velocity spectral slope: a v = —3.85 ± 0.02 


• The density spectral slope: a e = —3.3 (Stan- 
imirovic & Lazarian 2001) 
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